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Abstract 

We consider in detail the self-trapping of a soliton from a wave pulse 
that passes from a defocussing region into a focussing one in a spatially in- 
homogeneous nonlinear waveguide, described by a nonlinear Schrodinger 
equation in which the dispersion coefficient changes its sign from normal 
to anomalous. The model has direct applications to dispersion- decreasing 
nonlinear optical fibers, and to natural waveguides for internal waves in 
the ocean. It is found that, depending on the (conserved) energy and 
(nonconserved) "mass" of the initial pulse, four qualitatively different 
outcomes of the pulse transformation are possible: decay into radiation; 
self-trapping into a single soliton; formation of a breather; and formation 
of a pair of counterpropagating solitons. A corresponding chart is drawn 
on a parametric plane, which demonstrates some unexpected features. In 
particular, it is found that any kind of soliton(s) (including the breather 
and counterpropagating pair) eventually decays into pure radiation with 
the increase of the energy, the initial "mass" being kept constant. It is 
also noteworthy that a virtually direct transition from a single soliton 
into a pair of symmetric counterpropagating ones seems possible. An ex- 
planation for these features is proposed. In two cases when analytical 
approximations apply, viz., a simple perturbation theory for broad initial 
pulses, or the variational approximation for narrow ones, comparison with 
the direct simulations shows reasonable agreement. 
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1 Introduction 



Nonlinear spatially inhomogeneous waveguides give rise to a number of effects 
which are of interest by themselves, and also find important applications in 
such diverse fields as optical dispersion- decreasing fibers (DDF) [Q and natural 
waveguides for internal waves in the ocean with a shear-flow background Q . It 
is easy to understand that most nontrivial effects take place in the vicinity of a 
critical point, where the waveguide's dispersion or nonlinear coefficient changes 
sign. 

The critical points corresponding to the wave propagation of the nonlinear- 
Schrodinger (NLS) type were classified earlier in Ref. ||, where it was demon- 
strated that the most interesting one is that at which the sign of the dispersion 
coefficient a changes. However, unlike the case when the nonlinear coefficient 
in the corresponding NLS equation changes its sign, this case is not amenable 
to a consistent analytical consideration, hence systematic numerical simulations 
are necessary. 

In a very recent paper Q , the propagation of a wave pulse in this model was 
simulated for a situation when a changes its sign, in a self-focussing medium, 
from anomalous (admitting the existence of bright solitons) to normal (for which 
bright solitons do not exist). Accordingly, disintegration of an initial soliton-like 
pulse into radiation wave fields was considered. Despite the apparent simplicity 
of the process, a number of quite nontrivial features were found and qualitatively 
explained, the most interesting one being a double-humped structure in the 
region of the normal dispersion. 

For applications, particularly for those related to nonlinear optics of fibers 
and planar waveguides, especially relevant is the reverse process, i.e., formation 
of a soliton from a wave pulse crossing into the anomalous-dispersion region 
from the normal-dispersion one. Previously, this process was considered by one 
of the present authors in || in a purely analytical approximation, based on 
a variational technique. As we will demonstrate in this work, the system of 
ordinary differential equations derived in Ref. § from the underlying variable- 
coefficient NLS equation by means of the variational approximation, indeed 
provides for quite an accurate description of the pulse's dynamics in a parametric 
region where the approximation is relevant. Nevertheless, most results to be 
reported in the present work were produced by systematic direct simulations of 
the NLS equation. 

A modified NLS equation, valid near the zero-dispersion point for the nonuni- 
formly guided nonlinear wave propagation, was introduced in the previous works 

iu z + a(z)u u + 2\u\ 2 u — iSu ttu (1) 

where u is the local amplitude of the guided wave, z and t are the distance 
along the waveguide and the so-called reduced time (see, e.g., the derivation of 
the corresponding nonlinear Schrodinger equation for optical fibers in the book 
||), a(z) is the above-mentioned sign-changing variable dispersion coefficient, 
and S is the third-order-dispersion (TOD) coefficient which, generally, should 
be included in the case when the usual dispersion becomes very weak || . 

Here we consider solutions to Eq. (Ill) for a model with a continuous piecewise- 
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linear dispersion: 




a(z) = < z -K z < 1, (2) 



which takes into regard saturation of the dispersion after it has changed sign. 
This particular configuration can be easily realized in experiments with DDF 
and, generally, adequately represents the situation that we aim to consider. 

We will consider the evolution (for z > —1) of a pulse represented by the 
following natural initial configuration, 

u(z = -l,t) = Asech(ht). (3) 

Simulations demonstrate that the evolution of other smooth, unchirped initial 
pulses is very similar to that which is studied in detail below for the initial 
condition (||). 

2 Zero third-order dispersion 

If TOD is negligible, Eq. ([!]) takes a simpler form, 

iu z + ot(z)u t t + 2\u\ 2 u = 0. (4) 
Both equations (|l|) and conserve two quantities, viz., the "optical energy" 



and the momentum 



E= / \u(z,t)\ z dt, (5) 



P=- I {uu* t - u*u t )dt. (6) 
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These equations admit a Lagrangian representation. In particular, for the sim- 
plified equation (^) , the Lagrangian density is 

C = ^(u*u z ~ uu* z ) - a(z)\u t \ 2 + \u\ 4 . (7) 

Thus for zero TOD and a(z) chosen as per Eqs. (||) and (||), the pulse 
evolution is governed by the equations 

iu z + zu tt + 2\u\ 2 u — 0, (— 1 < z < 1) (8a) 

iu z + utt + 2\u\ 2 u = 0, (z>l), (8b) 

with 

u(-l,t) = Asech(ht). (8c) 

The evolution can be characterized by the two parameters, amplitude A and 
inverse size h of the initial pulse ([3|), or, alternatively, by its conserved energy 
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and initial "mass" . The "mass" , which is not a conserved quantity of the NLS 
equation, is defined as 

/oo 
\u(z,t)\dt. (9) 
-OO 

For the initial condition (g), M = M(z = -f) = nA/h, and E = 2A 2 /h. 

Since dispersion is constant for z > 1, the system (g) can be considered as 
an initial-value problem for the focussing NLS equation, with the initial condi- 
tion at z = 1. Thus the asymptotic solution will consist of a finite number of 
solitons represented by the discrete spectrum, and dispersive radiation repre- 
sented by the corresponding continuous spectrum. The number of solitons can 
be found via the Inverse Scattering Transform for the NLS equation by solving 
the Zakharov-Shabat (ZS) eigenvalue problem. To this aim, we define 

q(t)=u(z = l,t). (10) 

Assuming that \q\ decays as \t\ —> oo, the ZS eigenvalue problem is based on the 
linear equations 

iv t — Xv = qw, (H a ) 
iw t + Xw = qv (Hb) 

for auxiliary Jost functions u and v, the overbar designating the complex con- 
jugate. 

The spectrum of the eigenvalue A consists of its continuous part on the real 
axis and discrete eigenvalues, for which nontrivial solutions exist with functions 
v and w decaying exponentially as \t\ — > oo. The number of the discrete eigen- 
values in the upper half-plane is then equal to the number of solitons that will 
evolve from the wave packet q(t). These discrete eigenvalues, in general, have a 
nonzero imaginary part. 

Note that, for a symmetric q, which is the case here, if A is an eigenvalue, 
so also is —A. Solitary purely imaginary eigenvalues then correspond to a single 
zero-velocity NLS soliton, while multiple ones give rise to stationary breathers. 
If A has a nonzero real part, then the pair (A, —A) corresponds to a pair of 
counter-propagating solitons having velocities of equal magnitude but opposite 
sign. 

Thus, a two-stage process can be used to solve Eqs. (@). Firstly, Eq. ( |Sa|) 
is integrated from z = — 1 to z = 1, using a standard pseudo-spectral in t and 
fourth-order in z numerical method to obtain q(t) at z = 1. This is then used 
in Eqs. (|ll| ) to obtain the characteristics of the discrete spectrum from the 
ZS eigenvalue problem. Particularly, we want to know the total number of the 
discrete eigenvalues and the number of discrete eigenvalues on the imaginary 
axis. These two numbers can be found, using the methods elaborated in Refs. 
and |i. Introducing 

v — v cxp(— iXt), w — wexp(iXt), (12) 

we integrate Eqs. ( pd[ ) over a finite region t = [-L, L], where L>1, using the 
initial condition 0(— L) = 1 and finding the transmission coefficient v{L) = a(X). 
The discrete eigenvalues are then zeros of a(X) in the upper (Im{A} > 0) complex 
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half-plane. Since a is an analytical function of A, the number of the eigenvalues 
contained inside a closed curve T on the complex plane is 

To obtain the total number of eigenvalues, let T be the curve consisting of 
the real axis and a semi-circle with an infinitely large radius in the upper half- 
plane. On this infinite semi-circle, it can be shown that a(A) — > 1, thus the total 
number of eigenvalues is 

N =±r (i4) 

2m 7_ 00 o(A) V ' 

To obtain the number of discrete eigenvalues lying exactly on the imaginary 
axis, let r; m be the perimeter of a vanishingly thin but infinitely long rectangle 
enclosing the positive imaginary axis. Then the total number of the static (zero- 
velocity) solitons is N st = N(Ti m ). 

Before proceeding to full numerical solutions, we can consider two different 
analytical approximations for a solution of Eq. (|8a| ) in the interval z = [—1,1], 
valid, respectively, for long small-amplitude and short large-amplitude pulses. 
For the long pulses, we define the variables 

x = ht, u = Aip. (15) 

Then ([Sa| ) becomes 

i^p z + zh 2 ifj xx + 2A 2 \^\ 2 tP = 0, (16a) 

with 

ip(z = — 1, x) = sechx, (16b) 

Because we now want to consider a long small-amplitude pulse, we assume that 
h, A -c 1, such that 

h 2 = », A 2 = 7 V, (17) 

where 7 is an O(l) parameter and \i <C 1. We seek a perturbation solution in 
the form 

-0 = (Ro(x) + li 2 R 2 {z,x) + . . .)exp[ifji{if>i(z,x) + ^ 2 <t) 3 (z,x) + . . .)] . (18) 

Note that the 0(/x) corrections, Ri and cj>2, can be shown to be identically zero. 
The initial condition (16b) then gives 



At 0(h) we obtain 



R = sech(x). (19) 
Ro(j) lz = zR 0xx + 2 7 2 i?j], (20) 



and therefore 



ii = i(z 2 - 1)(1 - 2 sech 2 x) + 2 7 2 (z + 1) sech 2 x. (21) 
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Next, at 0(fi 2 ) it is found that R2 satisfies the equation 



R2z = ~z(2R 0x 4> lx + Ro<pixx), (22) 

and consequently 



i? 2 = 2 seclr x (4 - 5 sech x) 

\^ 1) - \{* 2 1) - 2 7 2 Q(^ 3 + 1) + ^ 2 - 1))) • (23) 

The amplitude obtained from this perturbation solution is compared with 
a numerical solution in Fig. |l[ As it is seen, the agreement is excellent. This 
perturbation solution would be expected to be valid for 

h 2 A 2 < 1 or E < 2(M a /n)i. (24) 

In the limit E — > 0, while Mq remains finite, it is apparent that at z = 1 the 
perturbation solution obtained by means of the above method is 

q(t) w Asech(ht), (25) 

in which case the variable-dispersion region is simply too short to have any effect 
on the very long pulse. For this initial condition, a well-known exact solution to 
the ZS eigenvalue problem was found by Satsuma and Yajima j^]. They showed 
that the solution consisted of N stationary solitons and a radiative component, 
where 



N = 



2Mn 



(26) 



[■] standing for the integer part. For Mq = Nn, where N is a positive integer, the 
radiation component is absent, and the solution consists purely of N interacting 
solitons. Examples of such solutions for Mq ^> ir can be found in Ref. |l(J. 

For short pulses, the variational approach proposed in Ref. || can be used. 
To this end, we assume an ansatz for the wave packet, 

u(z, t) = A(z) sech[</a(z)] exp [i(<f>(z) + b{z)t 2 )] , (27) 

where A, a, <f> and b are slowly varying real functions of z. Varying the averaged 
Lagrangian, 

/oo 
L dt, (28) 
-00 

with respect to these free parameters, it can be shown that A, <f) and b can be 
eliminated in favor of the pulse width, a, as 

^ = A( 2a A 2 ) = 0, (29a) 

az dz 

' ' (29b) 



Acta dz 
az oa z oa 
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Figure 1: A comparison of a numerical solution to Eq. (||) with the perturbation 
solution ( p"8| ) for Mo — 2n and E = 2, which correspond to p, = 1/16 and 
7 = 2 in Eqs. (15) and (17). The numerical solution and perturbation solutions 
are shown, respectively, by solid and dashed lines (that practically completely 
coincide), (a) The field \u\ at t = 0, and (b) |u| at z = 1. 



where a{z) is the variable dispersion coefficient in Eq. ( paj ) (a(z) = z in the 
case under consideration), and a(z) satisfies another variational equation, 



0. (30) 



dl\^dz~) + \n) V2^2~^3 

Since the initial chirp is zero, the appropriate initial conditions for (|30|) are 



, i da , 
a = h~ 1 , — = at z = -l. (31) 
dz 

In Fig. ||, an example of a numerical solution to Eqs. (|30|)-(^l|) is compared 
with a full numerical solution of (||) . The agreement between the two solutions 
is quite good, except for a narrow region near the maximum of \u(t)\ at z = 
1. Subsequently, the solutions experience rapid variations with z which the 
variational ansatz is unable to capture. 

Using these analytical approximations and direct numerical solutions, a few 
different types of the solutions to Eq. (||) can be identified, varying the initial 
mass and energy of the pulse. The results are summarized in Fig. |^. In the 
range considered, four particular types of the solutions were found. To construct 
this diagram, many numerical simulations were required, with particular atten- 
tion being taken to determine the boundaries between the different regions. Of 
course, these boundaries can only be determined to the accuracy shown here by 
using the two-stage numerical process described above. If only direct numeri- 
cal simulations are used, it becomes very difficult to gauge when a transition 
occurs. For example, a prohibitively large integration time is needed to dis- 
tinguish between the decaying-radiation and small-amplitude solitons, whereas 
this distinction can easily be made by considering the ZS eigenvalue problem. 
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Figure 2: A comparison of a numerical solution of Eq. (g) with the variational 
approximation (^) for Mq — tt/2 and E = 4. The numerical and variational 
solutions are shown by solid and dashed lines, respectively, (a) The field \u\ at 
t = 0, and (b) |u| at z = 1. 

For a small initial mass Mo, we find no discrete eigenvalues of the ZS scatter- 
ing problem, so the pulse completely decays into dispersive radiation, in which 
case | u | has a maximum at t = and is monotonically decaying with \t\. It 
is worth noting that near the R-C boundary, a double-humped structure can 
form at z > 1 with symmetric local maxima away from t = (which resembles a 
structure found recently in the same model with the reverse direction of the dis- 
persion change [9); however, eventually the pulse decays into a monotonically 
decreasing structure. An example of such a solution in shown in Fig. ^. For 
Mo — > 0, the evolution of the pulse at z > 1 is reasonably well approximated by 
the solution of the variational equation (|30|). 

As Mq is increased, two possibilities occur. In the first case, which occurs 
for small values of E, a single discrete eigenvalue appears on the imaginary 
axis. This corresponds to the "mass" of the pulse now being sufficiently large 
to form a soliton. Thus the asymptotic solution consists of a stationary soliton 
and dispersive radiation, see the example in Fig. g. 

The second case occurs for larger values of E. Here, for sufficiently large val- 
ues of the initial mass, the chirp has reduced sufficiently for the above-mentioned 
double- humped structure to develop into a pair of counter-propagating solitons. 
This corresponds to a pair of complex discrete eigenvalues appearing out of the 
real axis, where A = ±A r + iAj, A r and Aj being both real. An example is shown 
in Fig. |. 

In the last regime, a stationary breather, or a bound-state solution forms 
which asymptotically consists of a pair of interacting solitons. An example of 
this is shown in Fig. |t]. This can happen in two different ways. If the C — » B 
transition is considered, this occurs when the chirp has reduced sufficiently to no 
longer cause the splitting of the solitons. In terms of the ZS eigenvalue problem, 
this corresponds to the two eigenvalues approaching each other at a point on 
the imaginary axis. After their "collision" on the imaginary axis, two purely 
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Figure 3: A classification of numerical solutions to Eq. (g|) in terms of the 
energy E and initial mass Mo- R denotes the region of the radiative solutions, 
S the single-soliton solutions, B the bound-state (breather) solutions, and C 
the pair of counter- propagating solitons. 

imaginary eigenvalues appear, corresponding to two zero-velocity solitons. For 
the 5* — > B transition, this simply corresponds to the "mass" having increased 
sufficiently to give rise to a new zero- velocity soliton. This second soliton, like 
the one that forms in the R — > S transition, appears out of the origin and 
corresponds to a purely imaginary eigenvalue. 

Within the framework of the ZS eigenvalue problem, generic transitions that 
lead to the formation of a pair of counter-propagating solitons take place when 
two eigenvalues appear out of the continuous spectrum, such as occurs for the 
R — > C transition, or two imaginary eigenvalues collide and subsequently re- 
shape into a pair of complex eigenvalues, which occurs for the B — > C transi- 
tion. No evidence has previously been reported of a transition where a single 
imaginary eigenvalue splits into a pair of complex eigenvalues. Thus, from the- 
oretical considerations, we would expect a B region to always occur between 
the S and C ones, giving rise to a transition chain S — > B — > C. This inter- 
mediate B layer must disappear at a triple point, where the R and C regions 
become adjacent, as again the transition R — > C is admitted by the general 
analysis. For Al > 1.1 the S — > B — > C transition occurs as expected; however 
for M Q < 1.1 no evidence of the above-mentioned intermediate B layer could be 
found. Nevertheless, because a direct splitting of a soliton into a pair of sym- 
metric counterpropagating ones obviously contradicts the principle of continuity 
and therefore cannot take place in a model governed by a smooth differential 
equation, we conjecture that the missing intermediate B layer does exist, but it 
is too thin to be detected using the current numerical techniques. 
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Figure 5: The evolution of |w| for a single-soliton solution of Eq. (||) with E = 1 
and Mo = tt. Here, approximately 95% of the energy of the initial pulse is 
transmitted to the soliton. 
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Figure 6: The evolution of \u\ for a pair of counter-propagating solitons in Eq. 
(||) with E = 3 and Mo = 3tt/2. Here, approximately 93% of the energy of the 
initial pulse is transmitted to the two solitons. 
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Figure 7: The evolution of |u| for a bound-state (breather) solution of Eq. (||) 
with E = 2 and Mo = 37r/2. Here, between 90-95% of the energy of the initial 
pulse is transmitted to the breather. 
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If the initial "mass" of the pulse is fixed, say, to the value Mo = 1, it is 
apparent that as the energy is increased the number of solitons goes through the 
transition 1 — > 2 — > 0. This result appears counter-intuitive, as one would expect 
the number of solitons to monotonically increase with the energy. However, this 
transition can be explained by the influence of the pulse's chirp at z = 1. As 
E is increased, the strength of the chirp increases too. Thus, while for small 
values of E the chirp may have little effect on the pulse and one soliton is able 
to form, for larger values of the energy the chirp becomes strong enough to split 
the single soliton into a pair of solitons . For even larger values of the chirp, 
this pair is destroyed. A possibility could be to split the two solitons into four 
moving ones; however, because the mass at z = 1 is not sufficient to form four 
solitons, the pulse decays into radiation. 



3 Finite third-order dispersion 

In this section, we consider solutions to Eq. (|l|), supplemented by Eqs. (|^) and 
(|]), under the assumption that 5 is small (<5 <C 1) but finite. Therefore, the 
pulse evolution is now governed by the system 

iu z + zu t t + 2\u\ 2 u = iSuttt, (—1 < Z < 1) (32a) 
iu z + u tt + 2\u\ 2 u = iSuut, (z > 1) (32b) 



with 



u(z= -l,t) = Asech(ht). (32c) 

For z > 1, we are again dealing with a constant-coefficient equation, which 
is, however, no longer integrable (note that this equation, with both constant 
and periodically modulated coefficient in front of the usual dispersion term, has 
attracted a lot of attention in nonlinear fiber optics, see, e.g., Refs. Jjjl O, |l4| 
and references therein). Thus, only direct numerical simulations can be used to 
obtain the solutions. 

In the consideration of the transition from the anomalous to normal disper- 
sion (i.e., focussing to defocussing), which was the subject of Ref. the pulse 
width achieved a minimum at the zero-dispersion point z = 0, and consequently 
TOD was becoming more important as this point was approached. However, in 
the present case the pulse width achieves a maximum at z = and, in general, 
a minimum near 2 = 1. Thus, due to the the underlying assumption <5| <C 1, 
it would not be expected that TOD would be significant until well into the fo- 
cussing (anomalous-dispersion) region. Hence the effect of TOD on the pulse 
can be estimated by comparing the magnitude of the two constant-coefficient 
dispersion terms. If the pulse width is 0(a), then 

Suttt 



(33) 

From Figs. (|J)-(@) it is apparent that, for all the cases, a m ; n < ft, -1 , therefore 
TOD can be neglected if 

9 M 2 

\S\ « ||. (34) 
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Figure 8: The evolution of \u\ for the solution of Eq. @ with E = 2, M = tt/2 
and <5 = 0.01, demonstrating the effect of the third-order dispersion on the 
radiative solution which was shown, for the S — case, in Fig. ||. 

When TOD is taken into regard, it manifests itself in three ways. Firstly, 
TOD destroys the symmetry of the pulse. This is most dramatically seen if 
one considers the effect of TOD on the radiative solutions, an example of which 
is shown in Fig. || for a small value of the TOD coefficient, S = 0.01. In 
the defocussing (normal-dispersion) region (— 1 < z < 0), it is apparent that 
TOD has little effect. However, in the focussing (anomalous-dispersion) region 
(z > 0), it is apparent that TOD completely destroys the symmetry of the pulse 
well before the minimum of the pulse width is reached. Comparing this with 
Fig. [|, it is apparent that, nevertheless, the effect of TOD far away from the 
zero-dispersion point does not change the radiative type of the solution. 

The second effect of TOD, which is now well known, is that solitons are 
no longer localized, but rather generate a small-amplitude co-propagating os- 
cillatory tail (such nonlocal solitary waves have been termed "nanopterons" by 
Boyd Ijjl ) . Resonant generation of the tails in the context of the modified NLS 
equation including the TOD term has been considered in several works, see, 
e.g., Refs. |l(| [It]]. For a single soliton with speed c and amplitude A, it was 
demonstrated in Ref. Jl7j that, for 6 <C 1, the tail has a wavenumber 

k r = -6' 1 -c + 0{8), (35) 

and amplitude 

A^^exp^^), (36) 

where K w 8.58. For the single-soliton solutions considered in section 2, the tail 
generation is a major consequence of TOD (not shown here). 
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The symmetry breaking and appearance of the resonant oscillatory waves are 
clearly apparent if the effect of TOD on the counter-propagating soliton pair 
solutions is considered. An example is shown in Fig. ^. With S — 0, both soli- 
tons had the same amplitude and equal but opposite speeds. However, nonzero 
S destroys this symmetry. For positive 5, the leftward-propagating soliton is 
now reduced in amplitude and its speed decreases, whereas for the rightward- 
propagating soliton the amplitude and speed increase. The co-propagating os- 
cillatory waves (tail generation) can be seen both in the evolution of the wave 
field in Fig. ^|(a), and by examining its Fourier transform, which is shown (on 
the logarithmic scale) in Fig. ^(b) at z = 15. For the smaller soliton, the tail 
can be observed to the right of it, corresponding to the peak in the Fourier 
transform at k ss —3.5. For the larger soliton, it is more difficult to observe its 
tail in the evolution plot, however the front of this wave group can be seen start- 
ing to propagate rightwards from the pulse at z ss 2. These waves correspond 
to the well-pronounced peak in the Fourier transform at k rs —5.5 

Finally, we consider the effect of TOD on the breather that could exist at 
5 = 0. In this case, a major effect is splitting of the breather into two asymmetric 
solitons, provided that S is sufficiently large, an example of which is shown in 
Fig. no. Note that a possibility of splitting of a pulse under the action of TOD 



is also known in the constant-coefficient NLS equation, see, e.g., Ref. 12 . 

In the case shown in Fig. [l0[ the main breather-like soliton propagates to 
the right, but there is also a very weak wave propagating to the left. Whether 
the leftward-propagating wave is a soliton or just a packet of dispersive ra- 
diation cannot be determined. As in Fig. ||, there are oscillatory tail waves 
co-propagating with the main soliton. The front of these waves can be seen prop- 
agating rapidly to the right, and the peak in the Fourier transform at k —5.3 
(Fig. [H](b)) is associated with these waves. In simulations of the same case, 
but with 5 = 0.1 (twice as small, not shown here), the main breather soliton 
keeps a nearly zero velocity, and no small-amplitude leftward-propagating wave 
was detected. Qualitatively, this latter case, with the smaller value of 5, is very 
similar to what was shown in Fig. 

For each of the large-amplitude waves in Figs. (||) and (|Io|), a comparison 
of the observed co-propagating tail waves can be made with the predictions of 
Ref. fl7| , namely the expressions ( j35"[ ) and (j36|). This is shown in Table (|l|). 
For the small-amplitude soliton from Fig. |^, the comparison was not made, as 
the presence of the soliton with the larger amplitude and its tail renders the 
asymptotic analysis invalid for the smaller soliton. As can be seen from Table 
[j] for the dominant rightward propagating soliton, the agreement between the 
theoretically predicted and observed values of the amplitude is good; however, 
there is a conspicuous discrepancy, especially in the case shown in Fig. ^|, in the 
values of the wavenumber. Plausibly, the discrepancy can be fixed by including 
the next-order term in the expansion (]35|); however, one must also calculate the 
phase of the solitons in order to do this. 



4 Conclusion 

In this work, we have analyzed in detail self-trapping of a soliton from a wave 
packet that passes from a defocussing region into a focussing one in a spatially 
inhomogeneous nonlinear waveguide, described by a variable-dispersion NLS 
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Figure 9: The solution of Eq. @ with E = 3, M = 3tt/2 and 5 = 0.2, demon- 
strating the effect of the third-order dispersion on the counter-propagating 
soliton-pair solution, which was shown for case 5 — in Fig. ^. (a) The evolu- 
tion of the field \u\, (b) the absolute value of the Fourier transform at z = 15 
(note that the vertical scale in the panel (b) is logarithmic). 
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Figure 10: The solution of Eq. (||) with E = 2, M = 3vr/2 and (5 = 0.2, 
demonstrating the effect of the third-order dispersion on the bound-state solu- 
tion shown in Fig. (a) The evolution of the field |u|; (b) the absolute value 
of the Fourier transform at z = 23. 
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0.8 


2.3 


-7.3 


7 x 10~ 3 


-5.5 


6 x 10~ 3 




0.7 
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Table 1: A summary of predicted and observed amplitudes and wavenumbers 
for the oscillatory tail waves seen in Figs. (|) and @, where 5 = 0.2. A 
and c are, respectively, the soliton amplitudes and speeds, k r and A r are the 
wavenumbers and amplitudes predicted by Eqs. (35) and (|3^), while k a and A Q 
are the observed ones. 



16 



equation, in which the dispersion coefficient changes its sign from normal to 
anomalous. The model can be realized in terms of two (at least) very differ- 
ent but realistic physical applications: a dispersion-decreasing nonlinear optical 
fiber, and natural waveguides for internal waves in the ocean. It was found 
that, depending on the values of the (conserved) energy and (nonconserved) 
"mass" of the initial pulse, four qualitatively different outcomes of the pulse 
transformation are possible: decay into radiation; self-trapping into a single 
soliton; formation of a breather; and formation of a pair of counterpropagat- 
ing solitons. A chart of the corresponding parametric plane has been drawn, 
which demonstrates some unexpected features. One of them is that, with the 
increase of the energy, while the initial "mass" is kept constant, a soliton, a pair 
of the counterpropagating solitons, or a breather eventually decay into pure ra- 
diation. Another noteworthy feature is that a direct transition from a single 
soliton to a pair of symmetric counterpropagating ones seems virtually possible. 
An explanation for these features was proposed. In two cases when analytical 
approximations apply, viz., a straightforward perturbation theory for broad ini- 
tial pulses, or the variational approximation for narrow ones, comparison with 
the direct simulations shows a good agreement. 
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